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Abstract 

Changing the set of independent variables of Poincare gauge theory and considering, in a manner 
similar to the second order formalism of general relativity, the Riemannian part of the Lorentz 
connection as function of the tetrad field, we construct theories that do not contain second or higher 
order derivatives in the field variables, possess a full general relativity limit in the absence of spinning 
matter fields, and allow for propagating torsion fields in the general case, the spin density playing the 
role of the source current in a Yang-Mills type equation for the torsion. The equivalence of the second 
order and conventional first order formalism is established and the corresponding Noether identities 
are discussed. Finally, a concrete Lagrangian is constructed and by means of a Yasskin type ansatz, 
the field equations are reduced to a conventional Einstein-Proca system. Neglecting higher order 
terms in the spin tensor, approximate solutions describing the exterior of a spin polarized neutron 
star are presented and the possibility of the experimental detection of the torsion fields is briefly 
discussed. 

PACS: 04.50.+h; 04.20.Fy 



1 Introduction 

In general relativity, the procedure to consider, during the variation of the Lagrangian, the connection 
and the metric as independent of each other is known as Palatini, or first order, formalism. Indeed the 
classical Hilbert-Einstein Lagrangian contains no second derivatives of those fields. On the other hand, in 
the conventional metric approach of general relativity, where the only independent field is the metric, the 
Lagrangian contains second derivatives of the metric, although incidentally, the terms containing them 
can be eliminated by subtracting a boundary term. This is known as second order formalism. 

In Poincare gauge theory (PGT), the usual procedure is to consider the connection T ab m and the 
tetrad e° n as independent fields and consequently, the field equations are derived through variation with 
respect to those fields. For a detailed description of PGT, we refer to Refs. pQ and Usually, the 
Lagrangians are constructed from terms at most second order in curvature and torsion, and thus contain 
no second derivatives of the independent fields. We will therefore refer to this as the first order formalism. 

However, in view of the relation T ab m = T ab m + K ab m , where T ab m is the Christoffel connection, which 
can be expressed as function of the tetrad (and its derivatives) only, and K ab m is the contortion tensor, 
one might equally well consider K ab m and e^, as independent fields and carry out the variation with 
respect to those fields. Then, similar to the classical approach in general relativity (GR), the Christoffel 
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part of the connection is considered as a function of the tetrad. We will refer to this as the second order 
formalism. 

What could be the use of this formalism? Well, in conventional PGT, the most general Lagrangian 
is constructed under the requirement that it should not contain second and higher order derivatives of 
the independent fields and T ab m for the usual reasons (Cauchy problem). It is easy to see that, if 
we require instead the Lagrangian not to contain second and higher order derivatives of and K ab m , 
now seen as independent fields, then this will allow for other terms in the Lagrangian. Also, some terms 
previously allowed are now forbidden, as is the case, e.g., for terms of second order in the curvature, 
which will contain second derivatives of the tetrad field in this formalism. The Cauchy problem can 
now be stated with respect to and K ab m . Since the determination of the latter fields allows for the 
determination of T ab m too, no initial value problem will arise this way. 

On the other hand, one could argue that, in view of the underlying gauge structure of the theory, 
it appears quite unnatural to consider as fundamental variables a set different to the gauge fields of the 
Poincare group. It turns out, however, that both formalisms are completely equivalent. In other words, 
it is simply a matter of choice and convenience which formalism we use. This equivalence ultimately 
means that theories that contain higher derivatives in one formalism may be reformulated in the other 
formalism such that no higher derivatives occur. The change of the fundamental variables thus does not 
change anything concerning the viability or consistency of the theory, but it shows that certain theories 
that are apparently in conflict with the Cauchy initial value problem need not be so in reality. 

Summarizing, the second order formalism allows us to investigate different Lagrangians than those 
usually considered in PGT. This, in turn, will allow us to construct a theory with a complete classical 
general relativity limit, and thus a complete agreement with the experimental situation, without any 
constraints on eventual parameters of the theory, but with dynamical torsion fields arising in the presence 
of spinning matter. As we have pointed out in Ref. 0, this is not possible in conventional PGT for the 
following reasons. Since in the usual first order approach, the independent fields are the tetrad and the 
Lorentz connection T ab m , in order to get propagating modes for both fields, their first derivatives should 
appear quadratically in the Lagrangian. Therefore, apart from a term suitable of producing the general 
relativity limit, there has to be at least one term quadratic in the curvature tensor. In the presence of 
spinless matter, we wish the field equations to reduce to the Einstein equations of general relativity. If 
the Cartan equation (i.e., the connection equation) in this case leads to a vanishing torsion (Riemannian 
geometry), then the term quadratic in the curvature will necessarily contribute to the Einstein equation, 
which will thus not be of the GR form. (And therefore, in order to be in agreement with experiments, 
constraints will have to be imposed on the coupling constants of such theories.) The only exceptions to 
this are very artificial Lagrangians, with terms of the form R^R^ and R l uf. m -\Ri ^ km \ constructed from 
the antisymmetric part of the Ricci tensor and the completely antisymmetric (in the last three indices) 
part of the curvature tensor. Those terms will lead to contributions in the Einstein equation that vanish 
together with the torsion in the spinless case (in view of the Bianchi identities in the Riemannian limit). 
In the presence of spinning matter fields, those terms will eventually lead to propagating torsion, but in 
general, the field equations have very few and quite strange solutions, certainly not of the form expected 
for a Yang-Mills like theory. On the other hand, theories whose Cartan equation in the vanishing spin 
limit leads to a vanishing curvature (teleparallel geometry) and presenting a Yang- Mills like behavior for 
the Lorentz connection can easily be constructed using curvature squared terms, without problems with 
the general relativity limit. Those theories have been analyzed in detail in Ref. 0, where we have shown 
that they all present problems with an additional symmetry arising in the classical limit. Schematically, 



2 



we can write the linear approximation of such theories in the following form of wave equations 

DT ab rn = a ab m and Ug ik = T ik + 0{R), (1) 

which leads, in the spinless case, to T ab m = 0, and therefore the curvature corrections O(R) will vanish 
and we are left with the GR equation Dgi k — Ti k . The problem is that, even for T ab m fixed to zero, we 
can Lorentz rotate the tetrad field — » A a b e b n , without changing the metric gi k . Thus, the tetrad field 
is not uniquely fixed by the field equations, and therefore, the behavior of spinning test matter, e.g., the 
spin precession and the trajectory of a Dirac particle, cannot be uniquely predicted It is clear that 
this argumentation is valid independently of the formalism (first or second) one uses for the derivation 
of the field equations. It is a general problem of the set of equations r ab TO = and = T ik , i.e., of the 
teleparallel form of GR. 

On the other hand, a theory with an GR limit in a Riemannian geometry (vanishing torsion), should 
behave as 

UK ab m = a ab m and Ug ik = T ik + O(K), (2) 

where now the spinless case leads to a vanishing contortion K ab m — 0, and thus to a vanishing torsion, 
and therefore, as before, the second equation reduces to the field equations of GR. In contrast to (1), the 
set of equations (2) fixes the geometry completely. Indeed, if we know the contortion and the metric, 
than the independent field tensors, torsion and curvature, are fixed, up to a Poincare transformation of 
course, which is the underlying gauge symmetry of the theory [3] . 

However, as we have argued above, no Lagrangian of conventional PGT leads to equations of the form 
(2). In order to get such field equations using the conventional first order approach, one would have to 
use higher order derivatives of the fields (r m , e^J in the Lagrangian, which is (apparently) forbidden in 
view of the initial value problem. 

We will solve this problem using the second order approach and construct an explicit example of 
a theory with equations of type (2), containing only first order derivatives of the independent fields 
(K ab e a ) 

The article is structured as follows. In the next section, we will establish the equivalence of the first 
and second order formalism as far as the field equations are concerned. In section 3, we will compare 
the Noether currents and their conservation equations in the two approaches. Section 4 is devoted to 
the construction of an exemplifying theory that, using the second order formalism, allows for a full 
general relativity limit in the absence of spinning matter. Dynamical torsion fields will arise if the matter 
possesses a non- vanishing spin density. In section 5, extending slightly the Yasskin ansatz known from 
Einstein- Yang-Mills theory with internal symmetry group, we simplify the system of field equations for a 
special class of solutions, suitable for the description of spin polarized neutron stars. Finally, in section 
6, explicit solutions are discussed and in section 7, the effects of the torsion fields on the spin precession 
of a test body are briefly analyzed in view of an eventual experimental detection of the non- Riemannian 
contributions. 

2 Equivalence of first and second order formalisms 

Let us first give a short review of the basic concepts of Riemann-Cartan geometry and fix our notations 
and conventions. For a complete introduction into the subject, the reader may consult Refs. pQ and 0. 

Latin letters from the beginning of the alphabet (a, b, c. . .) run from to 3 and are (flat) tangent 
space indices. Especially, rj a b is the Minkowski metric diag(l, —1, —1, —1) in tangent space. Latin letters 
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from the middle of the alphabet k . . .) are indices in a curved spacctime with metric gik as before. 
We introduce the Poincare gauge fields, the tetrad e m and the connection T ab m (antisymmetric in ab), as 
well as the corresponding field strengths, the curvature and torsion tensors 

pab j^ab -pab , pa peb -pa -peb fn\ 

n lm — 1 rn,l 1 l,m ' 1 c/ 1 m 1 cm 1 ( l°J 

= e m,; — e ?,m + e m-T a b; — efr Q fcm . (4) 

The spacetime connection T l ml and the spacetime metric g^k can now be defined through 

a , -pa b api /r\ 

<e^77 ah = fflfc . (6) 

It is understood that there exists an inverse to the tetrad, such that e°e£ = S b . It can easily be shown 
that the connection splits into two parts, 

-pab pab , r^ab /tj\ 

1 m = 1 m + K m> ( 7 ) 

such that r ah TO is torsion- free and is essentially a function of e m . if ofc m is the contortion tensor (see 
below). Especially, the spacetime connection T l ml constructed from 

e m,l + F a bl e m — e i^ml (8) 

is just the (symmetric) Christoffel connection of general relativity, a function of the metric only. 

The gauge fields e m and T ab m are covector fields with respect to the spacetime index m. Under a 
local Lorentz transformation (more precisely, the Lorentz part of a Poincare transformation, see Ref. [3]) 
in tangent space, A a b (x m ), they transform as 

a \ a b pa a a a d r c A a a c /q\ 



The torsion and curvature are Lorentz tensors with respect to their tangent space indices as is easily 
shown. The contortion K ab m is also a Lorentz tensor and is related to the torsion through K l lm = 
U T l m + T m l - T ' Im). with = and analogously for T\ m . The inverse relation is T\ m = 

All quantities constructed from the torsion-free connection T ab m or Y\ will be denoted with a hat, 
for instance R l1 km = e l a e l b R ab km is the usual Riemann curvature tensor. 

We are now ready to compare the field equations arising in first and second order formalism and to 
establish their equivalence. 

Let C = Co + C m be the Lagrangian density in the first order formalism and C = Cq + C m the same 
Lagrangian in the second order formalism. The indices and m refer to the gravitational (or free) part 
and to the matter part respectively. Thus, we can write 

C(e m ,T ab m )=C(e ml K ab m ). (10) 

In concrete cases, the Lagrangians may differ by a total divergence, which does not influence the field 
equations and will therefore not be written out explicitely. For instance, in Ref. we considered the 
teleparallel Lagrangian 

e(R - -T m T lkl - X -T m T lkl + V ifc T™ ) = eR, (11) 
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where the left hand side is the first order form, while the right hand side is the usual Einstcin-Hilbcrt 
Lagrangian expressed in terms of the tetrad only. In (11), a divergence term has been suppressed. (In 
other words, the left and right hand sides are not really equal, but can be rendered equal by adding a 
physically irrelevant divergence.) 

The field equations in the conventional, first order formalism are well known, 

and -^=0. (12) 



As usual, they will be referred to as Cartan and Einstein equations, respectively. Both formalisms are 
related by equation (7), 

pab -pab la \ _i_ JS&b 

m 1 mV-"m) ' m' 

For the Cartan equation in the second order formalism, we find 

sc sc sr*, sc Q (13) 



6K ab m sr cd t 6K ab m sr ab 

Therefore, the Cartan equation has exactly the same form in both approaches. As to the Einstein 
equation, we easily find 

SC SC SC ST cd , SC SC ST cd n 

= 1 = 1 =0 f 14) 

Sef 5e™ ST cd t Se™ Se™ SY^ 5e™ ' v ; 

This differs by the second term from the corresponding equation in (12). This can be seen explicitely 
by carrying out the variation of this term using the explicit form of r cd ; from (8) and of the Christoffcl 
symbols t\ k . 

However, we can instead use immediately the Cartan equation SC/ST cd [ = (which we have shown 
to be equivalent to the Cartan equation in the second order formalism) to eliminate the second term in 
(14), and therefore we have the on-shell relation SC/Se™ = SC/5e™. 

To summarize, we have the following relations 

sc _ sc and sc sc 
s^~s^ an 6K« b m ~ sr ab m ' 

which shows clearly the equivalence of both sets of field equations. The above derivation also explains why 
some Lagrangians are allowed only in the first order formalism while others can only be considered in the 
second order formalism. The reason is that the terms resulting from the higher order derivatives (when 
the wrong formalism is applied to a certain Lagrangian) , can actually be eliminated using the second field 
equation, thereby bringing the field equations into a form that would have resulted immediately if the 
right formalism had been used right from the start. (By wrong and right, we simply mean the formalism 
where we have higher order derivatives, or we have not, respectively. Both formalisms are correct, since 
they are equivalent.) 

Note that in order to show the equivalence of both sets of field equations, we have to start with the 
full Lagrangian. We can of course not conclude that for given stress-energy tensor and spin density, the 
left hand sides of the field equations will have the same form in both procedures. The source tensors have 
to be modified at the same time. 
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Let us define the conventional canonical stress-energy and spin density tensors in the first order 
formalism 

T ^TeM and ° ah = e&r^> (16) 

as well as their second order analogues 

r^a I 5C m l SC m . . 

Ti= 2i74 and = e^- (17) 
Just as before (see (13) and (14)), one shows that 

1 sf cd 

°ab = CT ah and 1 i = 1 i + ^cd ■ ( 18 ) 

Now, use the relation (8), e° fc + Y a bk e\ = efT l ik , and the explicit form of the Christoffel symbols Y\ k to 
derive 

V cd — -\p ci p d — p di p c -4-p di p c — P ci p d 4-p di p ck P , p b —p ci p dk P , p b 1 

We can now evaluate T" in (18). The result is 

T~lm rrilrn , / _lmk kml klrn \ /in\ 
= J +2' _ _ ' * 

where ; denotes the usual covariant derivative with the Christoffel connection. 

Note that (19) looks basically like the Belinfante-Rosenfeldt relation between the symmetric Hilbert 
(or metric) stress-energy tensor T lm and the canonical (or Noether) tensor T lm known from purely 
Riemannian theory, see Refs. |J and [2], especially the discussion on relocalizations in Ref. 2 ). We will 
see that T lm is not yet the symmetric tensor, but it is a tensor that reduces to the symmetric Hilbert 
tensor for vanishing torsion. 

In order to prevent any kind of misunderstandings, let us emphasize once again that we do not 
modify the underlying gauge structure of the theory. Especially, the relation (10) makes clear the we do 
not change the coupling prescriptions in the matter Lagrangians. Thus, for instance, the Dirac particle 
still couples minimally to the gauge connection T ^ (which, from a mathematical point of view, is the 
fundamental field variable). The only difference is that we express this connection in terms of ef and 
K a \. Any change in the coupling prescription would lead to a fundamentally different theory, especially 
concerning the gauge symmetry. 

In the next section, we will derive the conservation equations for spin density and stress energy tensor 
in both approaches. 



3 The Noether identities 
3.1 First order formalism 

Let us briefly review the derivation of the Noether identities related to local, tangent space Lorentz trans- 
formations and to general coordinate transformations. Under an infinitesimal Lorentz transformation (9), 
with A a b = 5% + e a b (e ab = —s ba ) the independent fields e™ and r ab m undergo the following change: 

CT^afc ^ab pa r.b -nb ac i c m _ c m /or\\ 

° T m = - £ ,m ~ r cm £ ~ r cm e and 5e a = £ a e c • ( 2 0) 
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The first equation can be written in the short form ST ab m = —D m e ab . The change in the matter La- 
grangian therefore reads 

6C m = S -^6e™ + -^ST ab m = e(2Tl-] + D m0 — ) £(lc , (21) 



where we have omitted a total divergence. The covariant derivative operator D m is defined to act with 
T ab m on tangent space indices and with T l ki (torsion free) on spacetime indices. The requirement of 
Lorentz gauge invariance therefore leads to 

2T [ac] + Dm/ cm = Q ( 22 ) 

Slightly more complicated is the case of the coordinate invariance. Under an infinitesimal coordinate 
transformation, 

Z l =x l +£,\ (23) 

the fields T ab m and ef n transform as spacetime covectors (or one-forms), while the inverse tetrad e™ is a 
contravariant vector with respect to the index m. Thus, 

e™{x) = e™{x)+rA and f ab m (Z) = T ab m (x) - £ fc m r Qb fc . 

Since we are interested in the change of the Lagrangian under an active transformation, we have to 
evaluate the change of the fields at the same point x, and thus have to express the transformed fields in 
the old coordinates, i.e., 

e™(x) = e?(x) - e™ k (x)e, f ab m (x) = T ab m (x) - £*T ob mife (£). 

In the £-terms, we can replace e™(x) by e™(x) and T ab m (x) by r a6 m (x), since the difference will be of 
order £ 2 . Finally, we find 

8e? = CW-CW=G^"^ ( 24 ) 

sr ab rn = t ab m (x)-r ab m (x) = -e, m r ab k -er ab m , k . (25) 

Let us note at this point that, if we define the parameters e a = — e^£ m and e ab = £, k F ab k , the above 
transformations can be written in the form 

A P a — e a p c — n r a — J-T a 

Xrab t-v ab I nab 

01 m — U m E £ -ft lm- 

Since we are also free to perform an additional, independent, Lorentz transformation (20), we may again 
consider e a and e ab as independent of each other. In this form, the transformations are sometimes referred 
to as generalized Poincare gauge transformations, with e a generating the translations and e ab the Lorentz 
rotations pQ. We will remain in the usual interpretation of active coordinate transformations and return 
to the form (24)- (25). The change of the Lagrangian reads 

SC m = S J^6eJ a n + -^ST ab m 
5e™ a 5T ab 



= -2(eJ m e a),kt, -2eJ m e k £ +{ea ab 1 fe ), m £ -ecr ab 1 mM £, 



k 
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Requiring 5C m = and regrouping carefully the terms, we finally get 

(B m T m b )e b k + T m b T» mk = \R ab mk a ab m + ^T ab k (D m a ab m + 2T [ab] ). (26) 

The last term vanishes with (22) and the first two terms can be rewritten in a different form, using the 
relation T\ m = —2K\ m -, in the second term, and the explicit form of the covariant derivative D ln in the 
first term, to get 

rjim T^irnrp ^ T> a b ~r m f r > r 7\ 

1 k-m~ JX k 1 rni — -^ri mk a ab • \ z ' ) 

Recall that K %r \ is antisymmetric in im, so that classical matter (with Tj mi ] = and a ab m = 0) will 
satisfy the general relativistic conservation law T k ™ n = 0. 



3.2 Second order formalism 

In the second order formalism, we have to evaluate the changes of the fields e™ and K ab m under Lorentz 
and coordinate transformations. The Lorentz transformation now reads 

5e™ = e a c e™, 6K ab m = e a c K cb m + e b c K ac m , (28) 

which shows clearly that K ab n transforms as Lorentz tensor. Evaluating 
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and requiring SC m = we easily derive the following relation 

= 21 1 1 - K m a b +K m a b . (29) 

This does not look like a conservation equation for the spin density. We see that in the case where the 
torsion vanishes, the stress-energy tensor becomes symmetric. Therefore, we can see T m as a generaliza- 
tion of the symmetric Hilbert tensor. Let us note at this point that if we had considered as independent 
the fields K l \ and e a m (instead of K ab and e£j, we would have found a stress-energy tensor that is 
always symmetric. This, however, does not present any advantages (the variation of the gravitational 
Lagrangian would in most cases become more complicated) and so we preferred to stay a step closer to 
the conventional gauge approach, considering K ab m with its natural tangent space indices. The advantage 
of this approach is the linear relation between T ab m and K ab m and the resulting complete equivalence of 
the Cartan field equation in first and second order formalism. The alternative approach, with K l lm and 
e^, has been used in the framework of Einstein-Cartan-Dirac theory in the past @]. It leads of course to 
a symmetric Einstein equation. 

Under coordinate transformations, the fields transform exactly as in the previous section (because 
just as T ab m , K ab m too is a spacetime covector), i.e., 

6e? = C(*)-C(z)=f>* ( 3 °) 
6K ab m = K ab m (x)-K ab m (x) = ~e m K ab k ~eK ab k . (31) 
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In the same way as in the previous section, we derive the conservation equation in the form 



[y\ rpm\ a i rpmrpa m Jfli 

\ u mJ- a ) e kT 1 a 1 mk <J J« :m k 



1 

2 

-Y a K cb — Y b K ac -4- K a K cb A- K b K 

ck m ck m ' cm k ' cm k 



K ab — K ab 4- V a K cb 4- V b K ac 

k,m m,k ' cm k ' cm k 

^cb j^b r^ac , T^a r^cb , j^b r^ac 



We used (29) to put the equation into this form. A shorter way of writing this equation is the following 
(D m T m a )e£ + T™T a mk — —a li m . im K h k — - {R ab m k — R ab mk) <J ab n ■ (32) 
3.3 Discussion 

Comparing (22) and (29), we can derive the following relation between the antisymmetric parts of the 
stress-energy tensors: 

i2^rp[lm] rp[lrn]^ ^Imk (33) 

This relation is easily verified using the explicit form of T lm in terms of T lm , Eq. (19). 

Also, from (26) and (32), we get a relation between the divergences of the stress-energy tensors of the 
form 

{rjim rjim\ 1 r)db _ m /oa\ 

K 1 k - 1 k)\m - 2 ^ mfc°afc ■ (34) 

We used (33) to bring the relation into this form. Again, this relation can be derived directly from (19), 
using the symmetry properties of <j lmk and of the Riemann tensor R l lkm - 

Note that both sets of equations, (22,26) and (29,32), are valid, it is not the one or the other set. (The 
fact that the Lagrangian is invariant under Lorentz gauge and under general coordinate transformations 
does not depend on the choice of the independent field variables.) They are simply different equations 
for different quantities. The important thing is to use the right quantities as right hand side of the 
gravitational field equations (this depends on the choice of our formalism, first or second order) . 

Another point concerns the physical interpretation of T lk and T lk . If we take the point of view 
that the measurable quantities are not the densities, but rather the integrated quantities, namely the 
momentum vector and the spin tensor [5], then it is easy to see that T lk and T %k lead, in general, to 
different definitions of the momentum. Apart from this ambiguity, there is the problem that even from 
a single stress-energy tensor, one can define many inequivalcnt momentum vectors, integrating T"|, T l al 
T lk etc. (see Ref. [H] for an interesting discussion on this point). Usually, one would like to define 
the momentum in a way that is conserved. This is of course not possible in our case, since whatever 
momentum vector one defines, it will always be subject to gravitational forces coming from curvature 
and torsion, and therefore it cannot be conserved. (The other way around, whatever momentum vector 
one uses, it will always be conserved in the limit of vanishing curvature and torsion.) As it seems, it is 
rather a matter of convention how one wishes to define spin and momentum, under the only restriction 
that for vanishing fields, they should be conserved. 

In view of the rather unusual equation (29), which does not look like a conservation equation, one 
might be tempted to reject the set of current densities T a m , (J ah m and to claim that the canonical currents 
T a m ,a ab m as they arise in the conventional, first order formalism, should be considered as the true, 
physical currents. However, let us recall the fact that, for instance, the canonical stress-energy tensor 
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T a m of the Dirac Lagrangian contains non-axial torsion parts, which, as is well known, do not couple to 
the Dirac particle in the first place, and are neither contained in the Dirac equation, nor in the Lagrangian 
(where they cancel out with the hermitian conjugated term). We have brought this to attention in Ref. 
[5], to show that one has to be very careful on how to extract physical quantities from T a m . Thus, neither 
can the conventional, first order canonical tensors, be used without precautions to derive physically valid 
results. 

In order to gain at least some insight into the contents of T lk and T' tk , let us take a look at the flat 
limit of our equations, i.e., assume a spacetime with — rjn-. In this case, we find from (19) 

flva = T lm + _^lmk _ ^kml _ ^klmy ^ (35) 

If we define the momentum vector as we do in special relativity, 

P k = [ T\ dS h (36) 



then, for the alternative stress-energy tensor, we have 

h = J n dSt = Pk + lJ J - a l „: - a li m ), t dS,. (37) 

Clearly, the last term can be converted into a 4d volume integral over the divergence of the integrand, 
which vanishes however, due to the antisymmetry in il of the expression in brackets. Thus, we have 

h = Pk- (38) 

Of course, we also have (see (34) or (35)) 

1 k,m ~ 1 k,m- ) 

Usually, the relations (38) and (39) are said to show the equivalence of the two stress-energy tensors. 
(Same charges, same conservation law.) More generally, a change of the stress-energy tensor and the spin 
in the form 

rplm ^ rplm _|_ j^mlk (40) 

a km — * °km - { X km ~ X mk) (41) 

with X mkl = —X mlk is called relocalization of stress-energy and spin and is said not to affect physical 
quantities. 

It is easily seen that (35) is a relocalisation of T lm , without, however, the corresponding transformation 
of the spin density. Nevertheless, (26) and (29) coincide in flat space. The remaining relation is 



rprn rp 



k. 



2^ k\ a im ,1 ~ ^ il" 3 m — ^ ml" 3% I 

+ 2 (^"fc.m ~ K ab m,k + ^° t + ^"cm^"*)'' 'ab ' ■ 
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Unfortunately, we cannot express the equation for the spin precession (22) or (29) in a form indepen- 
dent of the stress-energy tensor. (Using (35) in (33) leads to an identity.) However, since the spin density 
is the same in both approaches, it will be subject to the same equations of motion, expressed in one or 
the other way. 

We therefore conclude that in the flat limit, gi k = rjik, the equations of motions for the current 
densities obtained in both pictures are identical as far as their physical contents is concerned. 

In curved spacetime, this is not the case anymore. Especially, we get two different equations for 
the stress-energy tensor and we have to decide, based on physical arguments, how to define define the 
momentum vector. 

Let us once again refer the reader to the extended discussion on canonical and metric stress-energy 
tensors and their relation via relocalizations in the general framework of metric affine theories in Hehl et 
al., 0. 

4 Model construction 

Comparing equations (27) and (32), we see that the following tensor 

— ( F? ab fi ab \ 

"'k mk) 

K ab — K ab A- T a K cb + T b K ac 

k.7n m,k * cm k ' cm k 

~ r a cfc A c6 , Tl — T b ck K ac m + K a cm K c \ + K b cm K ac k (42) 

plays, in the context of the second order formalism, a similar role as the curvature tensor R ah i m in the 
conventional approach. Essentially, F ab mk has the structure of a Lorentz Yang-Mills tensor, but with T ab m 
replaced by K ab m . The terms containing the Christoffel connection T ab m are needed to get the correct 
tensor behavior under a Lorentz transformation. 

Apart from F ab nk , which we will use to produce the dynamical, Yang-Mills like, behavior of the torsion 
field, we will need a term that guarantees the correct general relativity limit of the theory. In the second 
order formalism, this is simply the Einstein- Hilbcrt Lagrangian, expressed in terms of the tetrad field, R. 
For reasons that will become clear later, we also include a mass term K ab m K ab m into our Lagrangian. 

Thus, we start with 

Co = e(aR - bF ab mk F ab mk - cK ab m K ab m ). (43) 

This Lagrangian contains no second order derivatives of the independent fields ef and K a \, apart from 
the well known higher order terms contained in i?, which can, however, be eliminated by subtracting a 
surface term. On the other hand, if Co were to be treated in the first order formalism, there would arise 
higher order derivatives in the tetrad fields, coming from the term R ab i m R a b lm contained in the second 
term of (43). 

Including the matter Lagrangian C m , we are let to the following set of field equations 

- 46 D k F ab mk + 2cK ab m = a ab m (44) 

-aG kt = T M+ T W+T ( ll+f ki , (45) 

where 

t { H = -2b(F cb mt F cb m k -^g kl F cb ml F cb ml ) (46) 



ipab 
mk 



it 



r ( H = -c{K\K chl -\g kl K c \Kj) (47) 

t$ = -2b[(F lm £K l ™), + (Fj m K k n J,-(F k ^Kr n ), 

-{F l mkn K™ n h + (F km ?K lm n h - (F l m ?K k ™ n );i] . (48) 

We recognize in (44) a Yang-Mills type equation for the field K a b m , endowed with a mass term proportional 
to c. (Recall that we have defined D to act with the full Lorentz connection T ab m on tangent space indices, 
and with the Christoffel symbol on spacetime indices.) The corresponding symmetric, traceless stress- 
energy tensor is found in (46) and the stress-energy contribution from the mass term is given by (47). 

(31 

The only non-symmetric contributions in (45) are contained in the unusual term r y k ' i} a term that 
results from the F ab m couplings in F ab mk . In this term, the Lorentz gauge structure of the theory is not 
obvious, since the spacetime and Lorentz indices of F ab mk and K ab m are completely mixed up. This is also 
the reason why we chose to write down Eq. (45) immediately in terms of spacetime tensors, as opposed 
to the original form G a i = T" + . . . that results from the variation with respect to e* . In contrast to (44), 
which respects clearly the underlying gauge structure, nothing useful is revealed by keeping (45) in its 
mixed form, and it is convenient to pull down all quantities to physical spacetime such that one has to 
deal only with one kind of indices. 

Concretely, we find for the antisymmetric part of (45) 

f [ Jfc l = -2b{F^ m K kn m - F k n lm K m m ),i. (49) 

The calculations are rather lengthy, but one can directly verify that (49), together with (44), leads to the 
Noether identity (29). 

In the classical, spinless limit, it is immediately clear that (44) leads to the groundstate solution 
K ab m = 0, while (45) reduces to 

G% k — Ti k , 

which are the field equations of general relativity. (We took the opportunity to fix the parameter a to 
a = — 1 to be consistent with the conventional choice of units.) Also, in view of (19), there is no need, in 
the classical limit, to distinguish between T lk and T lk . 

Thus, the theory described by (43) allows for a complete general relativity limit in the absence of 
spinning matter, while a source with intrinsic spin, e.g., an elementary particle or a spin polarized 
neutron star, will give rise to dynamical torsion fields via (44). This cannot be achieved within the 
conventional first order formalism if one is not willing to allow for higher order derivatives, and therefore, 
no such theory has yet been presented. In order to avoid misunderstandings, let us emphasize that this 
does not mean that there are no consistent Poincare gauge theories based on the first order formalism 
without higher derivatives. What it means is that such theories do not have a full general relativity limit. 
Therefore, constraints on the free parameters of the theory will arise from the classical experiments, 
based on planetary motion and light propagation. This is rather disappointing, since it essentially means 
generalizing GR, but holding the changes small enough to ensure that they do not affect classical physics 
in a sensible way. On the other hand, in our model, classical measurements are not affected at all, 
independently of the coupling constants 6, c, and therefore, the post-general relativistic effects will only 
arise if spinning matter is introduced. Thus, the generalization affects only spinning matter, in view of 
which the theory was generalized in the first place. In short: Additional fields (torsion) for additional 
degrees of freedom (spin) . 

Finally, let us make a remark on the mass term in (43). The reason for its introduction is easily seen in 
context of the classical limit. Suppose c = in (44)-(45). Then, every GR solution, satisfying in addition 
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F ab mk = 0, will be solution to the field equations with cr Qb m = 0. However, from F ab mk = 0, we can not 
conclude to the unique solution K m — (we will see an explicit example in the next section). However, 
changing K ab m by a term that leaves F ab mk unaffected will have physical consequences (e.g., concerning 
the spin precession of a test particle) and therefore, this additional freedom cannot be accepted. This is 
not confined to the classical limit. From (43), it is obvious that, for c = 0, the field equations will only 
determine F ab mk , but not K ab m , allowing therefore, in general, for a whole class of solutions for K ab m . 
This will become very clear in the context of the solutions we will analyze in the next section. The 
introduction of the mass term obviously resolves this problem. 



5 Yasskin type reduction 

Since the classical limit of our theory coincides completely with general relativity, the only way to test 
the new features related to the presence of torsion fields is to look at field configurations generated by 
a source with intrinsic spin. One way is for instance to study interactions between elementary particles. 
On the astrophysical level, one can consider massive, spin polarized bodies, especially neutron stars, and 
analyze, e.g., the evolution of a system of two such bodies or, which is simpler from the computational 
point of view, study the motion and spin precession of a test body (a smaller neutron star or elementary 
particles) in the vicinity of a central source. 

In this section, we try to get solutions for the fields at the exterior of such a macroscopic spinning 
body. As discussed extensively in Ref . p] , such a matter distribution (located at the origin) is described 
by a spin density of the form 

C = W». (50) 
where p(x) is a suitably normalized density function and a a b is the integrated spin tensor. This spin 
density also coincides with that describing the generalized Weyssenhoff fluid in Riemann-Cartan space- 
times j7J. Then, following Yasskin jS], we look for solutions that have essentially the same structure as 
the current density, but in view of our rather special Lagrangian, we make the ansatz not for the Lorentz 
gauge potential T ab m , but for the contortion K ab m , i.e., we set 

K ab m = a ab A mi (51) 

where the field A m will have to be determined by the field equations. In Ref. with a similar ansatz, the 
Einstein- Yang- Mills field equations were reduced essentially to an Einstein-Maxwell system of equations, 
under the assumption that the gauge charges (in our case, <j a b) are constant. In order to preserve Lorentz 
invariance, one cannot require o a b to be constant, but instead has to use the constraint 

fW = 0, (52) 

which, with the help of (51), can also be written in the form 

D l( r Qb = 0. (53) 

In other words, the contortion fields, which in our second order approach play a role similar to the Yang- 
Mills potentials, drop out from the covariant derivative of a a b- This is similar to Ref. [S], where the 
gauge covariant derivatives coincide, for the specific ansatz, with the partial derivatives, which makes the 
requirement that the charges be constant a covariant one. 

What does Eq. (52) mean physically? Well, in view of a physical interpretation of the spin tensor, 
one will have to use further constraints, the so-called spin supplementary condition (SSC), like a^u 1 = 
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or similar. Roughly, there will be a local coordinate system where (52) reduces to — const and one 
might then choose the SSC <x;o = and define the spin vector c^ 1 = e^ x a u \ (u, A, v = 1, 2, 3). In other 
words, (52) is simply the covariant expression of the fact that the source is described by a constant spin 
vector. (Constant in the sense of a fixed vector, pointing in a fixed direction, as opposed to a radially 
orientated, spherically symmetric spin vector field, for instance.) This is just what we expect for a spin 
polarized body. Note, on the occasion, that a a bO' ab = 2a 2 . Therefore, let us define the spin s through 

a ab a ab = 2s 2 . (54) 

In practice, we will apply (52) only to the contortion field (51), and not to (50), since we will only 
deal with exterior solutions here. The above arguments are therefore only a simplified version of what 
should be the result of an analysis of the interior dynamics of the neutron star, which is beyond the scope 
of this article. For our purposes, the star can equally well be considered to be pointlike (i.e., p(x) = d(x) 
in (50)), as long as (52) is satisfied. 

With the ansatz (50)-(52), the field tensor (42) reduces to 

F ab k =a ab F ik) (55) 
where F- lk = A k ,i — Ai yk , whereas the field equations (44) and (45) take the simple form 

-4bF"% = -2cA rn + p(x)u m (56) 
G tk = -4bs 2 (F m F™-^g tk F lm F lm ) 

~2cs 2 (A k A l - ^g ik AiA l ) + f ik . (57) 

Thus, with a slight extension of the method presented in Ref. [jjj, we have brought our equations into an 
Einstein-Proca system. 

Several remarks are in order at this point. First, the ansatz (51) is supposed to hold at the exterior 
of the source. We do not know anything about the field configurations inside the source. Nevertheless, 
we have included the source terms in (56)-(57) in order to be able to fix eventual parameters of the 
solutions (like the Schwarzschild mass etc.). You may consider the source terms as boundary conditions, 
having in mind that the above form of the field equations is valid only at the exterior. Secondly, and not 
unrelated to the first remark, you may have noticed that equation (57) is symmetric, while one expects 
an asymmetric stress-energy tensor for a spinning body. This, of course, reflects again the fact that the 
equations do not hold inside the matter distribution. In order to ensure the continuity of the solutions 
at the boundary of the source, the stress-energy tensor should have a symmetric limit at that boundary. 
Recall however that the tensor Tik, as we have outlined at the end of section 2, is not the canonical stress- 
energy tensor, but a tensor that is already quite close to the symmetric Hilbert tensor. More precisely, 
we have shown that its antisymmetric part is given by (see (29)) 

21 = K m a b - K m a b ■ 

Using the spin density (50) and assuming that K ab m tends to the form (51) as we approach the boundary, 
we see that vanishes without any restrictions on the stress-energy tensor. Therefore, the symmetry 
of (57) is not a constraint on the stress-energy tensor, but simply the result of the specific form (51) of 
the torsion. 
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On the other hand, if we assume that K ab m has a similar form also at the interior of the source (this 
is to be understood as the mean value over a macroscopic region, not on an elementary particle level), 
then we can conclude that the stress-energy tensor is symmetric allover, and that it essentially coincides 
with the Hilbert tensor of general relativity. This is of course in contrast to conventional, first order, 
Poincare gauge theory, where a symmetric stress-energy tensor is only possible for vanishing spin. 

For the simplified equations (56)-(57), the reason for the necessity of the mass term (i.e., the c-terms) 
is especially clear. For c = 0, we are dealing with an Einstein-Maxwell system. Thus, the torsion will 
only be determined up to a gauge transformation A m — ► A m + / >m . However, measurable results, like 
the evolution of spin and momentum of a test particle in the given field configuration 5 , will depend on 
the full torsion tensor. Otherwise stated, the field equations for c = do not completely determine the 
physical fields. 

Unfortunately, the mass term, necessary on theoretical grounds, is rather disturbing from a compu- 
tational point of view, since the known solutions of the Einstein-Maxwell system do not easily generalize 
to the Proca field. 

6 Discussion and solutions of the field equations 

In the last section, we have reduced the field equations of our specific Poincare gauge theory to a purely 
Riemannian problem: The description of a massive vector field in curved spacetime. The discussion of 
the system (56)-(57) is therefore subject to the framework of general relativity, and can be found in 
the corresponding literature. Nevertheless, for the sake of completeness, we will briefly sketch the main 
features of this system and refer to the original articles for details. 

It is well known that in the cases of practical interest, no exact solution of the Einstein-Proca field 
equations has been found to date. Especially, no generalization of the spherically symmetric Reissner- 
Nordstr0m solution, nor of the axially symmetric Kerr-Newman solution are known for the case of a 
massive vector field. Moreover, it has been established [5] that there is no black hole solution other then 
for Ai = 0. This is a special case of the no-hair theorem for black holes. Here, we are interested in neutron 
stars. In view of the no-hair theorem, the vacuum solutions we eventually find for the exterior of the star 
will appear with a naked singularity at the origin. This, however, does not mean that they are unphysical. 
As vacuum solutions, they are restricted to the exterior of the star and will have to be completed by 
a suitable interior solution, thereby removing the singularity. On the other hand, the absence of black 
hole solutions indicates that there is a rather fundamental difference between the Maxwell and the Proca 
case. Especially, it seems that there are no solutions that tend to the Einstein-Maxwell black holes (Kerr- 
Newman or Reissner-Nordstr0m) for c — > (c is the mass parameter of the vector field in (56)-(57)), as 
might have been expected. 

A detailed analysis of the static, spherically symmetric case has been performed and numerical so- 
lutions have been derived in Refs. JOj-^^l- The authors of the three articles essentially agree in their 
conclusions, and there is no need to repeat the analysis here. The main feature of the solutions is that, for 
large distances from the source, the metric tends to the Schwarzschild solution. This is in agreement with 
the expected exponential decrease of the massive vector field known from the special relativistic limit. 
Indeed, for gik = riik, starting with the general spherically symmetric, static field Ai = ((p(r),i/j(r)-) in 
the comoving system of the source, u l = 5$, one easily shows from (56) that ip( r ) = 0j an d that tp(r) is 
solution to the equation 

- AbAip + 2ap = p, (58) 
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with solution 




In order to fix the proportionality factor, we have to take a look at equations (50) and (51). In contrast 
to the usual form of the Maxwell (or Proca) equations, the charge cr {, has essentially been removed 
from Ai in the ansatz for K ab m , i.e., in the flat case, the density p(r) is supposed to be normalized to 
/ p(r)d 3 x = 1. It is then easy to see (consider the special case c = and p(r) — 8(r)) that the complete 
solution reads 

^ = 16^ < 59 > 

Since the systematic numerical analysis can be found in Refs. |10 | -|12 | . we will apply here a different 
approach to the system (56)-(57). Namely, we will proceed in the spirit of the WKB approach of quantum 
mechanics, which consists in expanding the system in terms of Planck's constant. In our case, this means 
expanding in terms of the spin tensor a a b or simply in terms of the parameter s introduced in (54). To 
zeroth order (i.e., for s — 0), we simply get the Schwarzschild solution with Ai = 0. We are interested in 
the first order corrections to this solution. Thus, in (56)-(57), we simply neglect the terms in s 2 . Again, 
form (57), we find the Schwarzschild solution. (Just as in Einstein-Maxwell theory, where the corrections 
to the metric are of second order in the electric charge, here, the corrections are of second order in the 
spin.) At this stage, we are dealing with a Proca field on a Schwarzschild background. Finally, in order 
to solve equation (56), we expand the metric in terms of the Schwarzschild mass m and neglect the term 
of order (f(r)^ and higher, leaving us with the solution (59). 

Summarizing, at large distances, our solution looks like the Schwarzschild solution with Ai = 0. 
Approaching the source, the first order corrections are given by the solution (59) for Ai = (cp, 0, 0, 0). 
The next order of approximation consists in interaction terms of the gravitational potential m/r and the 
vector field Ai, leading to corrections in (59). Only at the next order (~ A 2 ) will the corrections to the 
Schwarzschild metric become apparent. 

Note that the Schwarzschild metric is not only a very good approximation because of the exponential 
decrease of (p, but also because of the absolute smallness of spin effects in general. In practice, rotational 
effects (leading to Kerr type corrections) will by far dominate the (intrinsic) spin effects. Not only will the 
rotational spin of a neutron star be much larger than its intrinsic spin, but moreover, the rotational effects 
lead to first order corrections in the metric, in contrast to the second order spin effects. Experimentally, 
this means that the spin effects cannot be observed by analyzing the geodesies, but rather by measuring 
directly the torsion field (through its interaction with spinning test bodies, see Ref. [5]). 

Let us emphasize that, although the Schwarzschild metric is subject to higher order corrections, it 
will retain its spherical symmetry to each order. The same holds for Ai, but not for K a \, which is given 
by (51). 

Unfortunately, the analysis of Refs. ^Hj-CSl has not been extended to the axially symmetric case. 
Since this is the case of practical interest (neutron stars are supposed to be subject to a rotation of very 
high frequency), we will extend the simple approach of the above considerations also to this case. As 
before, we thus neglect terms of second order in the spin. This, clearly, leads to the Kerr metric. Then, 
we also neglect interactions of the gravitational potential with the vector field Ai. 

First, let us recall the Kerr-Newman solution for the massless vector field, i.e., the axially symmetric 
solution of (56)-(57) for c = 0. Using spherical coordinates x l — (t,r,'&,<p), the Kerr-Newman metric 
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reads 



d S 2 = p2 - 2m I + q \ t> ^Ar 2 - p 2 dtf 2 
p 2 A 



(2mr — <7 2 )2asin 2 d , ,„ r9 9 , . , „ 2mr-o 2 . 1 9 , . 
^ didi? - [a 2 + r 2 + a 2 sin 2 -— - d<^ 2 , (60) 

with p 2 = r 2 + a 2 cos 2 A = r 2 — 2mr + q 2 + a 2 . Here, m, g, a are constants of integration. As is easily 
seen taking suitable limits, m is the Schwarzschild mass, a is related to the angular momentum of the 
source (see any textbook on general relativity) and the charge parameter q is proportional to the intrinsic 
spin s (see below). The corresponding field Ai has the form 

= W2b V a 2 cos 2 tf + r 2) ' 

The constant q can be evaluated by taking the Coulomb limit, i.e., take a — and compare with Eq. 
(59) for c = 0. The result is 

s\/2 
q ~ 16tiV&' 

and the solution takes the form 



(61) 



1 , r „ „ ar sin 2 1? , , . 

16-7TO p- 2 cos^ # + 

This is the exact solution for the massless vector field. A difference to the spherically symmetric case 
arises if we consider the flat limit of this solution. In the Reissner-Nordstr0m case, the solution Ai is 
simply the Coulomb potential, which is a solution of F % \ = in the flat case too. At first sight, it seems 
as if (62) is not a solution on the flat background ds 2 = At 2 — Ar 2 — r 2 (Ai) 2 + sin 2 dAtp 2 ). This, however, 
is the result of a false interpretation of the coordinate system. Indeed, the flat limit of (60) is found for 
q — > 0, m — > 0, i.e., 

ds 2 = d f _ r2+ f C °f" d Ar 2 _ ( r 2 + fl 2 cog 2 d)M 2 _ (fl 2 + ^ ^2 ^ 2 
r z _|_ a z 

which is easily shown to be flat (R l k i m = 0) and related to the Minkowski metric rfik — diag(l, — 1, —1, — 1) 
by the coordinate transformation 



x = \J r 2 + a 2 sin d cos ip 
y = V t 2 + a 2 cos $ cos ip 
z = rcosi). 

We now find that (62) is indeed a solution in the flat background expressed in the non-spherical coordinates 
r, i?, ip. Since it is rather difficult to generalize this to the massive case, we will neglect contributions that 
are of second order in the rotational momentum a. To that order, the coordinates used in (62) coincide 
with conventional spherical coordinates, and indeed the field 

1 1 1 „ „ a sin 2 i?. , . 

(-,0,0, ), (63) 



167T& r 
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which differs from (62) by terms of order a 2 and higher, is a solution to the Maxwell equations on a 
fiat background (in spherical coordinates). Physically, we see that (63) contains an electric monopolc 
and a magnetic dipole contribution, while (62) contains, as is well known, the complete series of odd 
electric multipoles (monopole, quadrupole. . . ) and even magnetic multipoles (dipole, octupolc. . . ). Thus, 
replacing (62) by (63), we essentially neglect quadrupole and higher moments. We then make the following 
ansatz for the massive vector field: 

1 , 1 „/ n „ „ a sin 2 # , ,, .„ 
^ = 16^ ( r /(r) '°'°' — 9{r)) - (64) 

Putting this into (56) yields, on a flat background, the unique asymptotically flat solution 

f(r) = e-^ ThT 

g(r) = (l + ^/2br)e-V^^ (65) 

where the constants of integration have been fixed by requiring the limit (63) for c — ► 0. 

Summarizing, the Kerr metric, as before the Schwarzschild metric, is a very good approximation for 
three reasons: (1) The spin is small compared to the angular momentum of the source. (2) Eventual 
corrections to the metric are of second order in Ai and thus in the spin (in contrast to the angular 
momentum contributions which are first order). (3) The field A i} and thus the metric corrections, 
decreases exponentially. 

The first order approximation of Ai is given by (64)- (65), neglecting interactions of the metric with 
Ai (order <~ Ai m/r) and second order terms in the spin and in the rotational momentum. 

Possibly, one can find exact solutions of (56) in Schwarzschild or Kerr backgrounds. This would yield 
the complete analogue of the WKB approximation, neglecting only terms of second and higher order in 
the spin, but including the spin-gravity couplings <~ Ai m/r. However, as we have argued before, it is 
not the scope of this article to analyze in detail the equations in the form (56)-(57), because this set of 
equations is well known from classical general relativity. 

Finally, one could also consider the possibility of a very small c, too small for the exponential decrease 
to be sensible on astrophysical scales. (In analogy to theories with a small photon mass.) In this case, 
one can consider the term in c as a gauge fixing term but irrelevant else. This means that we can simply 
take over the solutions from Einstein-Maxwell theory, i.e., the Reissner-Nordstr0m metric together with 
the Coulomb potential, or the Kerr-Newman metric with the field (62), but in the gauge that is consistent 
with c 7^ 0, which is essentially A\ — 0, as can be easily derived taking the divergence of (56). Note that 
(62) and also the Coulomb solution in its usual form are already in this gauge. Thus, the term in c forces 
us to choose the gauge A\ — 0, but the c-contributions will then be neglected in the field equations. 

Does every solution of (56)- (57) tend to a solution of the Einstein-Maxwell equations in the limit 
c = 0? Our approximate solutions seem to indicate that this is indeed the case. One should, however, be 
cautious, because there are nevertheless severe differences between the massive and the massless vector 
fields. One example is the previously mentioned absence of black hole solutions in the massive case. Thus, 
there might be surprises in the strong field region which have not been revealed by our approximations. 

Let us remind that the whole section was based on the specific ansatz (51), which was adapted to the 
macroscopic spin distribution (50). This does not mean that, even if we retain (50), there are no other 
solutions. Especially, the existence of black holes cannot be excluded a priori. On the other hand, for 
a different spin density, for instance the totally antisymmetric tensor of the Dirac particle, a different 
ansatz will be necessary, leading to entirely different solutions. 
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7 Spin precession in torsion field 



Finally, we address the question of the detectability of the torsion fields by analyzing the behavior of 
test bodies in a given field configuration. We confine ourselves to the concrete solutions of the previous 
section. 

It is well known that the equations for the spin precession and for the momentum evolution in a 
Riemann-Cartan geometry depend on the nature of the test body in question. Especially, elementary 
particles with different spin couple in a different manner to the torsion. A complete review, as well as a 
simple method for the derivation of the equations of motion, has been presented in Ref. 5 , where the 
references to the original articles can be found. 

For simplicity, we consider again the case where the test particle is a macroscopic, spin polarized 
body. Thus, the system (central source + test body) we are interested in is essentially a neutron star 
binary, under the restriction, for computational simplicity, that the central star has a much larger mass 
and can be considered to be at rest. 

The deviation from geodesic motion is given by a direct coupling of the test particles' spin to the 
curvature tensor. Even for the classical spin (i.e., the angular momentum in the body's rest frame), 
which presents a similar coupling to the curvature (Papapetrou equations in general relativity), these 
corrections have not yet been measured in experiment, and the contributions from the intrinsic spin and 
the non-Riemannian fields will be even smaller. Therefore, the only experimentally relevant equation is 
the spin precession equation, which reads |S] 

D^ 2) + f)S<£> = P t u k - P k u t , (66) 

where D denotes the covariant derivative with respect to proper time using the full connection Tj m and 
D the corresponding derivative with the Christoffel connection, i.e., Dc^' = dcr^/dr — T l im af^u m — 
rU4 2) « m and 6 4 2) = dS^/dr - f\ m S^u m - t{ m Sfu m . We denote by af k ] the (intrinsic) spin 

(2) 

tensor of the test body and by S^ k its angular momentum in the rest frame (due to rotation). The index 
(2) is attached to avoid confusion with the corresponding quantities of the source. Pi is the momentum 
vector of the test body, which is not necessarily parallel to the velocity u 1 . However, the right hand side 
of (66) is of higher order [S] and will be neglected in the following. Finally, we have to impose some 
conditions on the structure of the test body. It seems reasonable, for neutron stars, to make the ansatz 
(strong spin-rotation coupling) 

4 2) =5<4 2 \ (67) 

which leaves us with 

D<4 2) + ^K\ m <r$vr + J ±-K l km vg'>u m = 0, (68) 

where we have separated into Riemannian and non-Riemannian contributions. Next, introduce the spin 
vector ctLn = lrf kln ' 'ufcoj^ (jf klm = \g\~ 1 / 2 e lklm ) and use the spin supplementary condition of k u k = 
to find 

64) + ^—K\ m o\ 2) u m = 0, (69) 

where we have used the approximate validity of the geodesic equation Du 1 — and higher order terms 
have been omitted ■ This equation, together with the solutions of the previous section, can be used to 
evaluate the spin precession and especially to measure the influence of the torsion contributions. 
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Here, our scope is to illustrate the effects of the dynamical torsion fields. Therefore, we confine 

(2) 

ourselves to the simple case of the non-rotating test body (S^J = 0, i.e., g = 0) and use the spherically 
symmetric solution (no rotation of the source). This is not a very realistic case (as we have argued before, 
the rotational effects, when dealing with neutron stars, will in general be by far more important than the 
intrinsic spin effects), but it is a useful example to illustrate directly the post-general relativistic effects 
due to spin and torsion. 

With K a \ from (51), with A m = ((^,0,0,0) and introducing the spin vector of the source, o*,^ = 

i^%klm UkCr W^ usm g only the Newtonian order of the Schwarzschild background (compare with Ref. |S]), 
we find to lowest order for the relevant spin precession equation 

^ = l™[£x#V] + v {r)[dMx*% (70) 
at 2 r J 

where ip(r) is given by Eq. (59). The first term, with the orbital momentum L, is identical to the 
precession of the rotational spin in general relativity (as derived from the lowest approximation of the 
Papapetrou equations). This is a general result of Poincare gauge theory: In purely Riemannian space- 
times, no distinction can be made between intrinsic and classical (rotational) spin. The torsion, however, 
couples only to the intrinsic spin. Therefore, the other way around, only a particle with intrinsic spin 
can distinguish between Riemannian and non-Riemannian geometry. 

The above result is easily generalized to the more realistic case of rotating neutron stars. One can 
also proceed to a more complete treatment considering the neutron star binary as a two body problem. 
On the other hand, Eq. (70) is already enough to discuss the possibility of an experimental detection. 
It is clear that the spin effects are very small, even for completely spin polarized stars. Moreover, the 
torsion decreases exponentially, and practical results will strongly depend on the coupling constants b 
and c, especially on c. A possible, indirect, detection could be based on gravitational waves emitted by a 
merging binary system, since this will allow for conclusions on the evolution of the system in the strong 
field regime. 



8 Conclusions 

Considering, in the framework of Poincare gauge theory, the Christoffel part of the Lorentz connection 
as a function of the tetrad field in a formalism similar to the second order formalism used in general 
relativity and in metrical theories in general, we were able to present a theory that allows for dynam- 
ically propagating torsion fields, preserving nevertheless a full general relativity limit in the absence 
of spinning matter. This is not possible in the framework of the conventional, first order formalism, 
without introducing higher order derivatives of the independent field variables and apparently running 
into trouble with the Cauchy initial value problem. The equivalence of the second order formalism with 
the conventional approach has been established and therefore, ultimately, we have shown that in the 
conventional formalism, certain Lagrangians lead to consistent theories despite the appearance of higher 
derivatives. A concrete model has been presented, and for a specific class of solutions, the field equations 
have been reduced to an Einstein-Proca system using a slight modification of the Yasskin method known 
from conventional Einstein- Yang-Mills theories. Approximate solutions have been discussed for concrete 
cases, corresponding to the exterior of spin polarized neutron stars. In view of an eventual experimen- 
tal detection of the non-Riemannian effects, the spin precession equation of a test body in such a field 
configuration has been briefly analyzed. 
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